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Using data from RHIC we are able to systematically improve the two-component Dual Parton Model (DPM) 
event generator DPMJET-III. Introducing percolation parametrized as fusion of chains the model describes mul- 
tiplicities and pseudorapidity distributions in nucleus- nucleus collisions at all centralities. Guided by the d-Au 
data from RHIC we recalibrate the model to obtain collision scaling in h-A and d-A collisions. 



1. Introduction 

Hadronic collisions at high energies involve the 
production of particles with low transverse mo- 
menta, the so-called soft multiparticle produc- 
tion. The theoretical tools available at present 
are not sufficient to understand this feature 
from QCD alone and phenomenological models 
are typically applied in addition to perturbative 
QCD. The Dual Parton Model (DPM) [1] is such 
a phenomenological model and its fundamental 
ideas are presently the basis of many of the Monte 
Carlo implementations of soft interactions. 

One of the most prominent features of the Dual 
Parton Model has been so far the independent 
production and decay of hadronic strings. It 
seems, it is just this feature, which has to be 
modified in the extremely dense hadronic systems 
produced in central heavy ion collisions. 

2. Two component Dual Parton Model 

2.1. Hadron hadron collisions, the Monte 
Carlo Event Generator PHOJET 

PhojetI.12 [2,3] is a modern, DPM and per- 
turbative QCD based event generator describing 
hadron-hadron interactions and also hadronic in- 
teractions involving photons. Phojet replaces 

* Based on a poster submitted to the 17th International 
Conference on Ultra Relativistic Nucleus Nucleus Colli- 
sions, Jan. 11-17, Oakland, California USA 



the original Dtu jet model [4] , which was the first 
implementation of this combination of perturba- 
tive QCD and the DPM. 

The DPM combines predictions of the large 
N c , Nf expansion of QCD [5] and assumptions of 
duality [6] with Gribov's reggeon field theory [7]. 
Phojet, being used for the simulation of elemen- 
tary hadron-hadron, photon-hadron and photon- 
photon interactions with energies greater than 5 
GeV, implements the DPM as a two-component 
model using Reggeon theory for soft and lead- 
ing order perturbative QCD for hard interactions. 
Each Phojet collision includes multiple hard and 
multiple soft pomeron exchanges, as well as initial 
and final state radiation. In Phojet perturbative 
QCD interactions are referred to as hard Pomeron 
exchange. In addition to the model features as 
described in detail in [8], the version 1.12 incor- 
porates a model for high-mass diffraction dissoci- 
ation including multiple jet production and re- 
cursive insertions of enhanced pomeron graphs 
(triple-, loop- and double-pomeron graphs). 

High-mass diffraction dissociation is simulated 
as pomeron-hadron or pomeron-pomeron scatter- 
ing, including multiple soft and hard interactions 
[9] . To account for the nature of the pomeron be- 
ing a quasi-particle, the CKMT pomeron structu- 
re function [10] with a hard gluonic component is 
used. These considerations refer to pomeron ex- 
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change reactions with small pomeron-momentum 
transfer, \t\. For large \t\ the rapidity gap pro- 
duction (e.g. jet-gap-jet events) is implemented 
on the basis of the color evaporation model [11]. 
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Figure 1. Transverse momentum distributions of 
charged hadrons. The results of PHOJET (points) 
are compared to experimental data represented by 
lines, fitted to the data points. The data at y/s 
= 23 and 53 GeV are from the CERN-ISR[12], 
the data at 200 GeV are from the UAl Collabo- 
ration[13] and the data at 630 and 1800 GeV are 
from the CDF Collaboration[14]- 

The basic building blocks formed in these pro- 
cesses are color neutral strings. These strings are 
hadronized in Phojet using the Lund model as 
implemented in Pythia [15]. 

Phojet has been extensively tested against 
data in hadron-hadron collisions [8]. Further- 
more, in a number of papers LEP Collaborations 
compare many features of hadron production in 
7-7 collisions to Phojet, a rather good agree- 
ment is usually found. 

In Fig.l we compare charged hadrons according 
to Phojet against practically all data on trans- 
verse momentum distributions in p-p and p-p col- 
lisions from colliders. Please note, the points in 
this Figure are from the Phojet Monte Carlo, 
the data are represented by lines, fits to the data 
points. 



2.2. Collisions involving nuclei, the Monte 
Carlo Event Generator DPMJET-III 

The Dpmjet-III code system [16,17] is a Monte 
Carlo event generator implementing Gribov- 
Glaubcr theory for collisions involving nuclei, for 
all elementary collisions it uses the DPM as im- 
plemented in Phojet. Dpmjet-III is unique in 
its wide range of application simulating hadron- 
hadron, hadron-nucleus, nucleus-nucleus, photon- 
hadron, photon-photon and photon-nucleus inter- 
actions from a few GeV up to cosmic ray energies. 

The Gribov-Glauber Multiple Scattering For- 
malism: Since its first implementations [18,19] 
Dpmjet uses the Monte Carlo realization of the 
Gribov-Glauber multiple scattering formalism ac- 
cording to the algorithms of [20] and allows the 
calculation of total, elastic, quasi-elastic and pro- 
duction cross sections for any high-energy nu- 
clear collision. Parameters entering the hadron- 
nucleon scattering amplitude (total cross section 
and slope) are calculated within Phojet. 

Realistic nuclear densities and radii are used for 
light nuclei and Woods-Saxon densities otherwise. 

During the simulation of an inelastic colli- 
sion the above formalism samples the number 
of "wounded" nucleons, the impact parameter of 
the collision and the interaction configurations 
of the wounded nucleons. Individual hadron- 
nucleon interactions are then described by Pho- 
jet including multiple hard and soft pomeron ex- 
changes, initial and final state radiation as well as 
diffraction. 

As a new feature, Dpmjet-III allows the simu- 
lation of enhanced graph cuts in non-diffractive 
inelastic hadron-nucleus and nucleus- nucleus in- 
teractions. For example, in an event with two 
wounded nucleons, the first nucleon might take 
part in a non-diffractive interaction whereas the 
second one scatters diffractively producing only 
very few secondaries. Such graphs are predicted 
by the Gribov-Glauber theory of nuclear scatter- 
ing but are usually neglected. Further features 
of Dpmjet-III are a formation zone intranuclear 
cascade [21] and the implementation of certain 
baryon stopping diagrams [22]. 
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Dpmjet-III and earlier versions like Dp- 
MJET-II have been extensively tested against 
data in hadron-nucleus and nucleus-nucleus 
collisions[16-18,23]. The code is used for the sim- 
ulation of cosmic ray showers [24] . 

3. Comparing the original DPMJET III 
with RHIC data 

3.1. p p collisions at 200 GeV 

In Fig. 2 we compare the preliminary tt° trans- 
verse momentum distribution in p-p collisions at 
y/a = 200 GeV from PHENIX [25] to Phojet 
and find excellent agreement up to transverse mo- 
menta of about p± = 10 GeV/c. 
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Figure 2. Transverse momentum distributions 
of tv° -mesons produced in y/(s) = 200 GeV p- 
p collisions. The results of PHOJET are com- 
pared to preliminary experimental data from the 
PHENIX-Collaboration [25]. 

3.2. Au-Au collisions at 130 GeV in re- 
quiring Percolation 

The groups at Lisboa [26-28] and Santiago 
de Compostela [29] were the first to point out, 
that the multiplicities measured at RHIC are sig- 
nificantly lower than predicted by conventional 
multi-string models. A new process is needed to 
lower the multiplicity in situations with a very 



high density of produced hadrons like in central 
nucleus-nucleus collisions. The percolation pro- 
cess, which leads with increasing density to more 
and more fusion of strings, is one such mechanism 
[30,31]. 

Using the original Dpmjet-III with enhanced 
baryon stopping and a centrality of to 5 % we 
compare to some multiplicities measured in Au- 
Au collisions at RHIC for yfs NN = 130 GeV. 





N ch 


dN ch /dn\ri=o 


Dpmjet-III 


6031 


968 


BRAHMS [32] 


3860 ± 300 


553 ± 36 


PHOBOS [33] 




613 ± 24 


PHENIX [34] 




622 ± 41 



This comparison of Dpmjet-III with mul- 
tiplicity data from RHIC confirms: there is 
a new mechanism needed to reduce N c h and 
dN c h/dri\ ri= Q in situations with a produced very 
dense hadronic system. Therefore, in the next 
Section we will introduce percolation and chain 
fusion into Dpmjet-III. 

4. Percolation of hadronic strings in the 
modified DPMJET III 

We consider only the percolation with fusion of 
complete soft chains, these are chains where the 
transverse momenta of both chain ends is below 
a certain cut-off p^ slon . At the moment we use 
fusion _ 2 Q e v/c. The condition of percolation 
is, that the chains overlap in transverse space. We 
calculate the transverse distance of the chains L 
and K Rl-k and allow fusion of the chains for 
Rl-k < Rf usion . At the moment we use Rf^ on 
= 0.75 fm. This procedure is certainly the sim- 
plest approximation. A better procedure could be 
to introduce p± dependent transverse dimensions 
of the chains and then check for percolation. 

The chains in Dpmjet are fragmented using 
the Lund code Jetset as available inside the 
Pythia code [15]. Only the fragmentation of 
color triplet-antitriplet chains is available in Jet- 
set, however fusing two arbitrary chains could re- 
sult in chains with other colors. Therefore at the 
moment we select only chains for fusion, which 
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again result in triplet-antitriplet chains. Exam- 
ples for the fusion of two chains are: 

(i) A qi — q 2 plus a q 3 — (74 chain 

become a 51(73 — (72(74 chain. 

(ii) A q 1 - q 2 q 3 plus ag 4 -g 2 chain 

become a q\ 94 — q 3 chain. 
Examples for the fusion of three chain are: 

(iii) A q 3 - qiq 2 , a <j 4 - q\ plus a q 3 - q 5 chain _ 

become a (74 — (72*75 chain. 

(iv) A 94 — ^1, a q 5 — q 3 plus a q 5 — q 1 chain 

become a (74 — (73 chain. 
The expected results of these transformations 

are a decrease of the number of chains. Even if 

the fused chains have a higher energy than the 

original chains, the result will be a decrease of 

the hadron multiplicity Nhadrons- In reaction (i) 

we observe new diquark and anti-diquark chain 

ends. In the fragmentation of these chains we 

expect baryon-antibaryon production anywhere 

in the rapidity region of the collision. Therefore, 

(i) helps to shift the antibaryon to baryon ratio 

of the model into the direction as observed in the 

RHIC experiments. 

1200 



200 



DPMJET 0-5 % 
RHIC 0-5% 
DPMJET 5-10% 
RHIC 5-10% 
DPMJET 10-20% 
RHIC 10-20% 
DPMJET 20-30% 
RHIC 20-30% 
DPMJET 30-40% 
RHIC 30-40% 
DPMJET 40-50% 
RHIC 40-50% 





T I \ 




J:*' „ v 




















S "A. 




. . ■ h' 







Figure 3. Pseudorapidity distributions of charged 
hadrons in Au-Au collisions at y/s NN = 130 GeV 
for centralities 0-5 % up to 40-50 %. The points 
(with rather small error bars ) are from the Dp- 
mjet-III Monte Carlo with chain fusion as de- 
scribed in the text. With two exceptions data 
points are from the BRAHMS Collaboration [32]. 
The 77=0.0 points displayed at 77=0. 25 are from 
PHENIX [34] and the ?7=0.0 points displayed at 
77=-0.25 are from PHOBOS [33]. 



4.1. d-Au collisions in Dpmjet-III 

d-Au collisions are the first example where per- 
colation with fusion of chains is found to be im- 
portant. Pseudorapidity distribution of charged 
hadrons produced in minimum bias y/s — 200 
GeV d-Au collisions were measured at RHIC by 
the PHOBOS-Collaboration [35]. In Fig.4 we 
compare the PHOBOS data to Dpmjet-III cal- 
culations. Using Dpmjet-III without percolation 
with fusion of chains we find the Dpmjet distri- 
bution above the experimental data outside the 
systematic errors. Using Dpmjet-III with fusion 
of chains we find the Dpmjet distribution within 
the systematic errors. 
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Figure 4. Pseudorapidity distribution of charged 
hadrons produced in minimum bias ^fs — 200 
GeV d-Au collisions. The results of Dpmjet 
with percolation with fusion of chains are com- 
pared to experimental data from the PHOBOS- 
Collaboration [35]. We give at some pseudorapid- 
ity values the systematic errors as estimated by 
the experimental collaboration. 

Several RHIC experiments (see for instance 
[36]) find in d-Au collisions at large p± a nearly 
perfect collision scaling for ir° production. (Colli- 
sion scaling means Raa wl.O. ) The Raa ratios 
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are defined as follows: 

d 2 _jqA-A 



Raa = 



dp±drj 



jyA-A _ d 2 jyN-N 



(1) 



binary dp±dr) 



Here is the number of binary Glauber col- 

lisions in the nucleus-nucleus collision A-A. 

Dpmjet-III in its original form gave for ir° pro- 
duction in d+Au collisions strong deviations from 
collision scaling (Raa ~ 0.5 at large p±). The 
reason for this was in the iteration procedure to 
sample the multiple collisions in Dpmjet: some 
soft and hard collisions were rejected by this iter- 
ation procedure. Using a reordered iteration pro- 
cedure it was possible to obtain a nearly perfect 
collision scaling, see Fig. 5. 
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Figure 5. Raa ratio of ir° -mesons produced in yfs 
= 200 GeV d-Au collisions. The results of the 
modified Dpmjet are compared to experimental 
data from the PHENIX-Collaboration [36]. 

4.2. Au-Au collisions in Dpmjet-III 

Our Dpmjet-III studies on hadron production 
in central Au-Au collisions are still in a rather 
preliminary status. Dpmjet-III in its original 
form gave for 7r° production in central Au+Au 
collisions strong deviations from collision scaling. 
In contrast to the low p± behavior, the 7r° produc- 
tion at large p± was rather similar to the STAR- 
data[37]. Central to that behavior was that in 
the iteration procedure to implement the multi- 
ple collisions in Dpmjet-III some soft and hard 
collisions were rejected. 



The RHIC data on hadron production in d-Au 
collisions made it clear, that this iteration pro- 
cedure had to be changed in such a way that no 
essential deviations from collision scaling occur in 
d-Au collisions. The same corrections to the it- 
eration procedure then also influence Au-Au col- 
lisions in a way, that the deviations from collision 
scaling are reduced. Hence we have to recalcu- 
late all hadron production at large p± in Au-Au 
collisions and subsequently we have to introduce 
nuclear modifications to parton distributions and 
interactions of the scattering partons in the dense 
medium. These modifications have not yet been 
implemented so far. 

The only thing we can present at the mo- 
ment is a recalculation of large p± ir° produc- 
tion with Dpmjet-III with full percolation and 
chain fusion and with the kinematical corrections 
to change the iteration procedure. In Fig. 6 we 
present such Dpmjet-III results and compare 
them to the RHIC data[37]. As expected, at 
larger p± Dpmjet-III produces a p± distribution 
above the data. 

5. Elliptic flow in Dpmjet-III 

The elliptic flow v 2 is defined as a Fourier co- 
efficient of the single particle distribution [38,39] 

v 2 =< cos2((j) - $_r) > . 

is the orientation of the reaction plane and 
the angular brackets denote an average over many 
particles belonging to some phase-space region 
and over many collisions having approximately 
the same impact parameter. The elliptic flow v 2 
is a sensitive probe of the dense matter produced 
in high-energy nucleus-nucleus collisions. Data 
from RHIC on v 2 were published by most RHIC 
Collaborations [40-43]. 

Quite a number of different methods have been 
proposed to extract the elliptic flow from data 
or from Monte Carlo models, it is not the pur- 
pose of this contribution to give a overview over 
these methods. Let us only mention two meth- 
ods, which we use to analyze elliptic flow from 
Dpmjet-III. One method uses in its simplest 
form two-particle cumulants [44] . The second me- 
thod determines flow using Lee- Yang zeroes [45] . 
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Figure 6. Transverse momentum distribution of 
7T° -mesons produced in central yfs = 200 GeV 
Au-Au collisions. The results of the modified 
DPMJET-iU are compared to experimental data 
from the PHENIX-Collaboration [37]. 



We concentrate on Au-Au collisions at V s / 
= 130 GeV, for which experimental data for the 
integrated flow and for p± dependent flow were 
presented in Ref . [43] . In a multichain model like 
Dpmjet-III with non-interacting hadronic strings 
fragmenting into hadrons and for large secondary 
multiplicities, we expect of course a vanishing el- 
liptic flow. But there are two effects in Dpm- 
jet-III which are expected to change this picture. 
The first effect is connected to the percolation of 
hadronic strings in Dpmjet-III. An interaction of 
hadronic strings can be expected to also introduce 
a force between strings yielding some genuine el- 
liptic flow. The second effect are jets and minijets 
in Dpmjet-III. The number of large p± jets and 
minijets at RHIC energies is not large, the jets 
produce significant azimuthal fluctuations, which 
leads to non- vanishing w 2 especially at large p± . 

In Table 1 we collect for two centralities (and 
the pseudorapidity acceptance of the experiment) 
some integrated flow values calculated with both 



Table 1 

Integrated flow: v-i results obtained from Dp- 
mjet-III using two different methods (I: using 
two particle cumulants[44]; II: using Lee- Yang ze- 
roes [45]; III: gives the spurious flow according to 
method II) to STAR data [43] from RHIC. 



Ccntrality 


particles 


I 


II 


III 


STAR 


53-77 % 


charged 


0.0107 


0.0140 


0.047 


0.07 


53-77 % 


7r± 


0.0182 


0.0135 


0.051 




53-77 % 


if± 


0.060 


0.037 


0.164 




0-5 % 


charged 


0.00685 


0.00644 


0.0129 


0.02 


0-5 % 




0.00595 


0.00693 


0.0139 




0-5 % 




0.0188 


0.0215 


0.0419 





methods mentioned above. In Fig. 7 we com- 
pare the pj_ dependent elliptic flow obtained 
from Dpmjet-III using the two-particle cumu- 
lant method [44] for two different centralities to 
STAR data [43] from RHIC. 

Let us discuss the results of the comparisons in 
Table 1 and Fig.7. 

(i) The integrated flow values calculated from 
Dpmjet-III with methods I and II agree reason- 
ably well, however the spurious flow (column III 
of Tab. 1) calculated for method II indicates, that 
method II using Lee- Yang zeroes is not able to 
calculate the flow reliably. Also we find, that we 
are not able to calculate a significant p± depen- 
dent flow using method II. 

(ii) The integrated flow values in Table 1 calcu- 
lated from Dpmjet-III are significantly smaller 
than the experimental data from STAR, also we 
see in Fig.7 the p± dependent flow at low p± val- 
ues (which dominate the particle production) is 
significantly smaller than the experimental data 
from STAR. We conclude: the flow introduced in 
the model by the percolation and fusion of chains 
is not enough to explain the data. 

(iii) We find in Fig.7 that Dpmjet-III at large 
p± values shows larger flow than the data, this 
flow results from the jets and minijets as dis- 
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cussed above. This flow component in the experi- 
mental data is significantly smaller, we know that 
large p± jets in the Au-Au data are strongly sup- 
pressed, the present version of Dpmjet-III does 
not contain any mechanism for jet quenching. 
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Figure 7. Elliptic flow V2 as function of p±_ in 130 
GeV Au-Au collisions at RHIC. We compare for 
charged hadrons the Dpmjet-III results for two 
centralities to experimental data from the STAR- 
Collaboration [43] 



6. Summary 

The data obtained at RHIC are extremely use- 
ful to improve hadron production models like 
Dpmjet-III. 

Of particular importance are data on hadron 
production in p-p collisions, d-Au collisions and 
peripheral Au-Au collisions, in all of these colli- 
sions (unlike central Au-Au collisions) we do not 
expect any change in the reaction mechanism, 
which could not be accommodated into the mech- 
anisms as implemented in Dpmjet-III. 



Indeed, comparing Dpmjet-III to RHIC data 
we find two important corrections to be ap- 
plied to Dpmjet-III, which otherwise do not 
completely change the independent chain frag- 
mentation model: (i) Percolation and fusion of 
chains, the data from RHIC allow to determine 
the amount of percolation to be implemented 
into Dpmjet-III. (ii) Collision scaling of large p± 
hadron production in d-Au collisions: The data 
indicate that we have to change the iteration pro- 
cedure (of the selection of all soft and hard chains 
in nuclear collisions) in such a way, that collision 
scaling is obtained. 
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